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Abstract 



This paper describes a search for the Standard Model Higgs boson using data from e"'"e~ 
colhsions collected at center-of-mass energies of 161, 170 and 172 GeV by the OPAL detector at 
LEP. The data collected at these energies correspond to integrated luminosities of 10.0, 1.0 and 
9.4 pb~^, respectively. The search is sensitive to the main final states from the process in which 
the Higgs boson is produced in association with a fermion anti-fermion pair, namely four jets, 
two jets with missing energy, and two jets produced together with a pair of electron, muon or tau 
leptons. One candidate event is observed, in agreement with the Standard Model background 
expectation. In combination with previous OPAL searches at center-of-mass energies close to 
the resonance and the revised previous OPAL searches at 161 GeV , we derive a lower limit 
of 69.4 GeV for the mass of the Standard Model Higgs boson at the 95% confidence level. 
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1 Introduction 



Locally gauge-invariant theories of the electroweak interaction introduce spontaneous symmetry 
breaking to allow some of the gauge bosons to acquire mass while keeping the theory renormal- 
izable. The Standard Model (SM) is the simplest such theory and uses the self-interaction 
of a single doublet of complex scalar fields to produce spontaneous symmetry breaking. 
This model predicts the existence of one physical scalar particle, the Higgs boson, H°, whose 
couplings are fixed but whose mass is not predicted. 

Despite a wide experimental effort, the Higgs boson has not yet been discovered. The 
experimental lower limits for its mass, mno, obtained from large samples of TP boson decays, 
are published in 0. Recently the OPAL Collaboration has updated its result by including the 
data collected during summer 1996 at a center-of-mass energy of -\/s=161 GeV yielding a lower 
limit of 68.5 GeV on the Higgs boson mass @]. During the autumn of 1996 the center-of-mass 
energy of the LEP e+e~ collider was upgraded to 172 GeV. 

At these center-of-mass energies, the main production process for the SM Higgs boson is 
e"'"e~— >-Z*^H°. The dominant decay is H°— i>bb, with a branching ratio of approximately 86%. 
Other relevant decay modes are: H°— i>r"'"r~ (8%), H°— >cc (4%), and H°— i>gluons (2%) 0. For 
Higgs boson masses of current interest, these branching ratios exhibit only a mild dependence 
on the Higgs boson mass. 

The searches described address the principal final state topologies, which account for about 
95% of all Higgs boson final states, namely: (i) the four jets channel, e"'"e~— >Z'^H'^— >-qqbb; (ii) 
the missing energy channel, mainly from e"'"e~— >Z'^H'^^z/z/bb, with a small contribution from the 
W"''W~ fusion process e^e~ ^vv'R^^ (iii) the tau channels, e^e^^Z°H°^r+r^qq and qqr+r^; 
(iv) the muon and electron channels, predominantly from e"'"e~— >-Z°H°— ^-z^^/x^qq and e"*"e~qq, 
with the latter including a small contribution from the TPTP fusion process e"'"e~— >-e"'"e~H''. 

The present paper describes the analysis of the data collected at 170 and 172 GeV energies, 
an improved analysis of the four jets channel for the 161 GeV data based on likelihood rather 
than cut-based methods, and the derivation of a new mass limit. 



2 Detector, Data, and Simulations 

This analysis uses 1.0 pb^^ of data recorded with the OPAL detector at -^5=170 GeV, 
9.4 pb"^ at Vs=172 GeV and 10.0 pb"^ at ^=161 GeV for the four jets channel. 

OPAL is a multipurpose apparatus with nearly complete solid angle coverage and excellent 
hermeticity. The central tracking detector consists of two layers of silicon microstrip detectors 
with polar angleQ coverage |cos^| < 0.9, immediately outside the beam-pipe, followed by 
a high-precision vertex drift chamber, a large-volume jet chamber, and z-chambers, all in a 

^OPAL uses a right-handed coordinate system where the -\-z direction is along the electron beam and where 
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uniform 0.435 T axial magnetic field. A lead-glass electromagnetic calorimeter is located out- 
side the magnet coil, which, in combination with the forward calorimeter, gamma catcher, and 
silicon-tungsten luminometer provide acceptance down to 24 mrad from the beam direction. 
The silicon-tungsten luminometer serves to measure the integrated luminosity using small-angle 
Bhabha scattering events |^ . The magnet return yoke is instrumented with streamer tubes for 
hadron calorimetry and is surrounded by several layers of muon chambers. Events are recon- 
structed from charged-particle tracks and energy deposits ("clusters") in the electromagnetic 
and hadronic calorimeters. The tracks and clusters are required to pass a set of quality re- 
quirements similar to those used in previous Higgs boson searches . In calculating the total 
visible energy and momentum, ii^vis and Pvis, of events and of individual jets, corrections are 
applied that reduce the effect of double-counting of energy in the case of tracks and clusters 
associated to them. 



The signal detection efficiencies and accepted background cross sections are estimated using 
a variety of Monte Carlo samples all processed through a full simulation of the OPAL 



detector. The HZHA generator |T^, including initial-state radiation effects, is used to simulate 



Higgs boson production processes. The generated partons are hadronized using JETSET [|Tj]. 
Signal samples are produced for fixed values of mno between 40 GeV and 80 GeV. The estimates 
of the different background processes are based primarily on the following event generators: 
PYTHIA (Z77*->qq(7)), grc4f ^ (four-fermion processes including W+W" and Z°Z*), 
BHWIDE ll (e+e-(7)), KORALZ (/i+/x-(7) and t+t-{j)), PYTHIA, PHOJET ^ and 
Vermaseren (hadronic and leptonic two-photon processes). 



3 The Four Jets Channel 



The process e"'"e~^Z''H''— >qqbb amounts to approximately 60% of the Higgs boson signal 
topologies. It is characterized by four energetic hadronic jets, large visible energy and sig- 
nals from b-hadron decays. The backgrounds are Z°/7*-^qq with and without initial state 
radiation accompanied by hard gluon emission as well as four-fermion processes, in particular 
e"'"e~— >-W"'"W~. The suppression of these backgrounds relies on the kinematic reconstruction 
of the Z° boson and on the identification of b-quarks from the Higgs boson decay. The tagging 
of particles containing b quarks proceeds by detecting displaced secondary vertices in three 
dimensions exploiting the high-resolution obtained from OPAL's silicon microvertex detector 
and using leptons with high transverse momenta with respect to the jets to which they are 
assigned. 

The selection of candidate events and the suppression of background is done in two steps. 
A preselection using cuts is first performed in order to retain only events which have some 
similarity to the signal. The remaining events are then analyzed using a likelihood technique. 

This channel was already analyzed and published using cut-based methods for the 161 GeV 

+x points to the center of the LEP ring. The polar angle, 9, is defined with respect to the +z direction and 
the azimuthal angle, 0, with respect to the horizontal, +x direction. 
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data [Q], however, using likelihood methods and a better b-tag improves the signal efficiency 
significantly with the same expected background. The analysis of this channel was therefore 
repeated for the 161 GeV data with the new likelihood method. 

The analysis proceeds in the following way: 

First, large parts of the Z°/7*— i>qq background are eliminated by selecting well defined four 
jet topologies using the cuts listed below: 



^1) The events must qualify as a hadronic final state according to ref. ||20 



(2) The radiative process e"'"e — >Z''7^qq7 is largely eliminated by requiring that the effective 
center-of-mass energy, obtained by discarding the radiative photon from the event 
], is at least 140 GeV (150 GeV) at = 161 GeV (v^ = 170 - 172 GeV). 



(3) The final state particles of an event are grouped into four jets using the Durham algo- 
rithm |]22|. The jet resolution parameter, 1/34, at which the number of jets changes from 
3 to 4, is required to be larger than 0.005. 

(4) The Z°/7*— ^qq background is further suppressed by requiring that the event shape pa- 
rameter C ^3l, which is large for spherical events, is larger than 0.45. 



(5) All four jets are required to contain at least two tracks and two electromagnetic calorimeter 
clusters. 

(6) To discriminate against poorly reconstructed events, a kinematic fit, using energy and 
momentum conservation constraints (4C-fit), is required to converge with a probability 
larger than 0.01. 

The e"'"e~^H''Z° hypothesis is tested by a kinematic fit which, in addition to the energy 
and momentum conservation constraints, also forces two of the four jets to have an invari- 
ant mass equal to the Z° boson mass (5C-fit). This fit is applied in turn to all six possible 
associations of the four jets to the Z° and H° bosons. The fit is required to converge for 
at least one combination with a probability of at least 0.01. The combination yielding 
the highest x^-probability is selected. 



Table |I| shows the number of events selected in data and Monte Carlo at each stage of 
the cuts for both center-of-mass energies. After removing most of the two-photon background 
(cut(l)), a good agreement between the observed data and the expected background can be 
seen. 

Next, a likelihood technique is employed in order to classify the remaining events as either 
Z°/7*— i>qq (1), a four-fermion process (2) or Z°H°-^qqbb (3). To select signal events with 
low background, the information from 14 quantities (described below) which provide a good 
separation between the three different event classes is combined. Half of the variables are 
related to the kinematics of the events, and the other half are related to b-tagging. For each 
event the measured values Xi (i=1...14) of these variables are compared with probability density 
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Cut 


Data 


Total bkg. 


qq(7) 


4-ferm. 


Efficiency (%) 




161 GeV 








= 65 GeV 


(1) 


1572 


1389.4 


1345.9 


52 


99.9 


(2) 


395 


377.5 


351.1 


26.4 


89.5 


(3) 


65 


54.1 


38.0 


16.1 


81.4 


(4) 


51 


40.6 


26.2 


14.4 


75.6 


(5) 


49 


33.2 


21.5 


11.7 


70.5 


(6) 


26 


24.2 


14.0 


10.2 


62.3 


> 0.9 





0.75 ± 0.08 


0.49 


0.26 


32.1 



Cut 


Data 


Total bkg. 


qq(7) 


4-ferm. 


Efficiency (%) 




170-172 GeV 








mno = 68 GeV 


(1) 


1409 


1306.7 


1189.7 


117.0 


99.7 


(2) 


367 


381.2 


312.8 


68.4 


87.8 


(3) 


93 


84.7 


33.5 


51.2 


79.3 


(4) 


77 


70.9 


22.9 


48.0 


75.8 


(5) 


68 


60.3 


18.1 


42.2 


70.2 


(6) 


56 


50.3 


12.6 


37.7 


61.4 


C"^ > 0.925 


1 


0.88 ± 0.07 


0.34 


0.55 


28.4 



Table 1: The number of events after each cut of the selection for the data at ^/s — 161 (170- 
172) GeV and the expected background in the four jets channel. The background estimates 
are normalized to 10.0 (10.4 ph~^). The last column shows the selection efficiencies for the 
Z°H°— >qqbb final state. 
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functions normalized to unity obtained from Monte Carlo events processed through the full 
detector simulation. The likelihood for each event to belong to any of the three event classes 
is calculated as follows. For a single variable, the probability for an event to belong to class j 
is given by 



fiixi) 



where fHxi) denotes the probability density for event class j and variable i. The likelihood 
function for class j is defined as the normalized product of the pj(xj) 



J2k=i nr=i Pii^i) 

where n = 14 is the number of variables. In order to select an event as a candidate event its 
likelihood for being a signal event is required to be larger than a certain value depending on 
the preferred signal/background ratio. 

The first set of variables entering the likelihood function exploit the different kinematics of 
the background and signal events: the smallest angle between any pair of jets (1), the smallest 
dijet mass (2), the highest and lowest jet energy (3-4) and the dijet mass closest to mzo (5). 
All these quantities are calculated after the 4C-fit. The probability of the best kinematic 
fit requiring energy and momentum conservation and equal dijet masses (6) and the larger 
(3 = p/ E factor of the two dijet momenta calculated from the jet pair combination closest to 
the WW hypothesis (7) are also used. The best jet pair combination under the WW hypothesis 
is determined by minimizing (pdijoti — 2;)^ + (Pdijot2 — a^)^, where x is 30 GeV at 170-172 GeV and 
6 GeV at 161 GeV center of mass energy. This quantity vanishes for perfectly reconstructed 
on-shell WW events. 

The second set of variables is used to tag b-fiavored hadrons. Secondary vertices are iden- 
tified in each jet separately using the three dimensional extension of the method described in 



2^ . To improve the quality of the vertices, each accepted vertex is required to have at least 
two tracks containing 2 silicon microvertex hits in both r0 and z. The variables to discriminate 
between b flavor and lighter flavors are: sum of the decay length significances of the vertices 
found in the candidate Higgs jets (8)0, the maximum number of tracks with significantly large 
impact parameter (9)0, the invariant mass of the tracks in the vertex with the largest decay 
length significance (10), the smaller of the two Higgs jet masses (11), the sum of the momenta 
of the highest momentum track associated with each vertex of the two Higgs jets (12), the 
largest transverse momentum of an identified lepton (electron or muon) with respect to the 
corresponding jet axis (13) and the sum of the two largest numbers of significant tracks in the 
four jets (14). 

Finally the likelihood for each event is required to be larger than 0.90 at ^/s = 161 GeV 
and larger than 0.925 at = 170 - 172 GeV (i.e. C^{x) > 0.90 (0.925)). The likelihood 

^The decay length significance is defined as L/aL, wliere L is tlie tfiree dimensional distance between the 
primary and the secondary vertex and is the corresponding error. 

•^Significant tracks are defined by d/ad > 2.5 where d denotes the two dimensional impact parameter of the 
track and ad the error on the impact parameter. 
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requirement was tightened for the higher center of mass energy in order to retain the same level 
of expected background for both energies. 

The analysis was tuned with a reference mass of mno = 65 GeV (68 GeV ) at y/s= 161 GeV 
(170-172 GeV) with a resulting efficiency of 32.1% (28.5%). The resulting expected background 
is 0.49 (0.34) events from Z/'j* and 0.26 (0.55) events from 4-fermion processes. Other sources 
of background are negligible. The total expected background is 0.75±0.08(stat)±0.25(syst) 
(0.88±0.07±0.18) events. The background systematics are dominated by the error on the 
modeling of the variables used in the likelihood. 

For all fourteen input quantities a good agreement between OPAL data and the Monte 
Carlo (MC) distributions is observed. This can be seen in Figure |1] where some distributions 
of data and simulated background as well as a simulated 68 GeV Higgs signal are shown. The 
agreement between MC and data, in particular for b-tag related quantities, was also checked 
at y/s=mzo with data taken for calibration immediately before the high energy runs. 

Figures || (a) and (b) show the distributions of the signal likelihood for the preselected 
events from the 161 GeV and 170-172 GeV data. The shaded area shows the expectation 
from 4-fermion events, the white area is Z/'y* events and the dotted line represents the total 
background with a signal from a 65 GeV Higgs at ^/s = 161 GeV and a 68 GeV Higgs at y/s = 
170-172 GeV added. It can be seen that the expected signal is concentrated at large values 
of the likelihood. The peak of the background at large likelihood values is due to irreducible 
signal like four-fermions and QCD four jets events containing b flavor. 

No event survived the likelihood cut of 0.9 at ^/s = 161 GeV. With a cut at 0.925 in the 
likelihood variable one candidate event is selected from the data collected at ^/s = 172 GeV. 
The likelihood of this event is 0.993. The invariant mass of the two jets associated with the 
Higgs boson decay is 75.6 ±3.0 GeV. Both of the Higgs candidate jets have displaced secondary 
vertices with significances of 6.1 and 2.0, and charged multiplicities of 4 and 3. One of the jets 
assigned to the Z°-candidate has also a displaced secondary vertex with a significance of 4.6 
and multiplicity of 2. The event is shown in Figure 

The likelihood method has been also checked using calibration events at a LEP energy of 
^/s = rrizo where the physics is well established. Events classified as four jets are selected and the 
tracks momenta and clusters energies are scaled up to 172 GeV center of mass energy. The b-tag 
information and most of the kinematic event properties are preserved by the procedure. The 
likelihood values of the scaled events show a good agreement between data and MC simulation 
(see Figure ^ ). 

The signal detection efficiencies at ^/s = 161 GeV (170-172 GeV) are affected by the follow- 
ing uncertainties: Monte Carlo statistical error 4% (2%); description of preselection variables 
5% (2%); modeling of the variables used in the likelihood selection, 10% (7%). The latter was 
dominated by b-tag variables uncertainties which were estimated by varying the tracking reso- 
lution. The modeling was also checked using a re-weighting procedure. The difference between 
data and MC distributions was minimized by assigning weights to all MC events in an iterative 
procedure and recalculating the likelihood with the weighted MC distributions. The binning of 
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reference distributions yielded a contribution to the systematic error of 1% (2%). Taking these 
uncertainties as independent and adding them in quadrature results in a total systematic error 
on the signal efficiency of 12% (8%) relative error. 

Two additional analyses were performed as a check of the likelihood method. The first was 
a cut-based method similar to the published analysis of the 161 GeV data 0. This analysis 
gave, for similar residual background, signal efficiencies which were lower by 10-20% than those 
of the current analysis. Applied to the 170-172 GeV data, the cut-based analysis did not select 
the candidate event selected by the current selection because it failed the cut in the event shape 
parameter C. The C-parameter of this event is 0.546, close to the cut value of 0.55. 

Another multidimensional analysis utilizing an Artificial Neural Network (ANN) with 9 
b-quark tagging variables and 12 kinematic variables was also performed. Using the same 
Monte Carlo signal and background samples described previously, the ANN analysis achieved 
similar efficiency, signal-to-noise and systematic errors, to the likelihood approach. The ANN 
identified one Higgs boson candidate with a cut on the ANN output variable corresponding to 
a background of 0.9 events. This event was the same as the one identified in the main analysis. 



4 The Missing Energy Channel 



The e^e"— s>z/z/H°^z/z/bb process amounts to approximately 18% of the Higgs boson signal 
topologies with a small contribution (0.7% for mH0=68 GeV) from the W"^W~ fusion pro- 
cess e^e~ ^vvYl^^vuhh. These events are characterized by large missing momentum and 
two energetic, acoplanar, hadronic jets containing b-hadrons. The dominant backgrounds are 
mismeasured Z°/7*^qq events, four-fermion processes with neutrinos in the final state, such 
as W+W~^£^z/qq and W^e^i/— *-qq e^v with the charged lepton escaping detection and, in 
general, events in which particles go undetected down the beam pipe such as e^e~— >Z°7 and 
two-photon events. For most of these backgrounds, the missing momentum vector points close 
to the beam direction, while signal events tend to have missing momentum in the transverse 
plane. 

The event selection proceeds as follows: 



^1) To reduce two-photon and beam-wall interactions, the number of tracks passing the qual- 



ity cuts |jTO[ should be greater than six, and should exceed 20% of the total number of 
tracks in the event. The energy deposited in the forward detector, gamma catcher and 
silicon-tungsten luminometer must be less than 2 GeV, 5 GeV, and 5 GeV, respectively. 
The fraction of energy deposited in the region | cos 6*1 > 0.9 must not exceed 50% of the 
total visible energy in the event. The total transverse momentum of the event, P^g, must 
be greater than 1 GeV and the visible mass must satisfy m^is > 4 GeV. 

(2) To remove backgrounds in which particles go undetected down the beam pipe, the polar 
angle, 6*111188, of the missing momentum (Pmiss = — -Pvis) must satisfy | cos6'iniss| < 0.9. The 
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z component of the visible momentum, P^^g, is required to be less than 30 GeV. 

(3) The remaining two-photon background is eliminated by requiring > 8 GeV. As a 
precaution against large fluctuations in the measured hadronic energy, P^^ is recalculated 
excluding hadronic calorimeter clusters and is also required to be larger than 8 GeV. 

(4) The remaining events are reconstructed as two-jet events using the Durham algorithm. 
The axes of both jets are required to have a polar angle satisfying | cos 6*1 < 0.9, to ensure 
good containment. 

(5) The remaining background from Z°/7*— i>qq is characterized by two jets that tend to be 
back-to-back with small acoplanarity angles (where the acoplanarity angle is defined as 
180° — 0jj where is the angle between the two jets in the plane perpendicular to 
the beam axis), in contrast to signal events in which the jets are expected to have some 
acoplanarity angle due to the boost of the Higgs boson. This background is suppressed 
by requiring that the jet-jet acoplanarity angle be larger than 5°. 

(6) Since the Higgs boson recoils against a Z° boson decaying into a pair of neutrinos, the 
signal has a missing mass close to mzo. The remaining backgrounds, predominantly 
from well contained multi-hadron and four-fermion events including semi-leptonic W"*'W~ 
decays, typically have small missing masses. These backgrounds are reduced by the 
missing mass requirement 76^ < m^jgg < 120^ GeV^. The distribution of the missing 
mass squared after cut (5) is shown in Figure H. 

(7) W"''W~ events with one of the W bosons decaying leptonically and the other decaying 
into hadronic jets are rejected by requiring that the events have no isolated leptons. 

In this context, leptons are low-multiplicity jets with one (two or three) tracks associated 
to electromagnetic or hadronic energy clusters, confined to a cone of 5° (7°) half-angle, 
having an invariant mass less than 2.5 GeV and momentum in excess of 5 GeV. The 
lepton is considered isolated if the sum of the track energy and the electromagnetic energy 
contained between the above lepton cone and an isolation cone of 25° half-angle does not 
exceed 10% (15%) of the lepton energy. If the lepton cone has only one track, the isolation 
cone is not allowed to contain another track. 

(8) The remaining background is mostly from semi-leptonic W+W^ and Wez/ events where 
the charged lepton goes undetected down the beam pipe. These events are suppressed 
by requiring b-hadrons in the hadronic jets of events with a visible mass close to the W 
boson mass. Each of the two hadronic jets is required to contain a secondary vertex with 
at least two tracks, each containing two silicon microvertex detector r0 hits. The two 
vertices are reconstructed according to the three dimensional extension of the method 
used in P3. The correlation between the sum of the decay length significances, S^, and 



the visible mass of the event is shown in Figure ^ where the two-dimensional cut employed 
is also shown. 

The numbers of observed and expected events after each selection cut are given in Table ^ 
along with the detection efficiency for a 68 GeV Higgs boson. No events survive the cuts while 



11 



the expected background is of 0.55±0.05(stat)±0.10(syst) events. The background systematics 
is dominated by the error on the modeUing of the cut variables. 

One event is rejected only by the b-tag requirement. This event has two well-contained 
hadronic jets and evidence for an isolated muon pointing at the very forward direction where 
no tracking information is available. The muon appears as electromagnetic and hadronic clus- 
ters aligned with two muon chambers hits. With a well-measured muon, this event would be 
classified as W"*"W~— >qq/ii/ and rejected by the isolated leptons veto. 

The detection efficiencies as a function of the Higgs boson mass are listed in Table ^ 
These include a small correction (2.3%) due to accelerator-related backgrounds in the forward 
detectors which are not simulated. The detection efficiencies are affected by the following 
uncertainties: Monte Carlo statistics, 2.2%; modeling of the cut variables other than b-tagging, 
3.8%; b-tagging and uncertainties from fragmentation and hadronization 3.0% and lepton tag, 
5.6%. Taking these uncertainties as independent and adding them in quadrature results in a 
total relative systematic uncertainty of 7.7%. 



Cut 


Data 


Total bkg. 


qq(7) 


4-ferm. 


77 


Efficiency (%) 
mno = 68 GeV 


(1) 


2984 


2829 


713 


105 


2011 


83.3 


(2) 


1468 


1486 


302.1 


86.7 


1097 


73.6 


(3) 


173 


177.5 


121.0 


56.0 


0.37 


71.9 


(4) 


163 


165.9 


113.5 


52.0 


0.30 


64.3 


(5) 


53 


58.5 


18.3 


40.0 


0.30 


62.3 


(6) 


2 


2.2 


0.6 


1.5 


0.1 


55.2 


(7) 


1 


1.6 


0.52 


1.0 


0.1 


52.8 


(8) 





0.55±0.05 


0.275 


0.275 





42.7 



Table 2: The numbers of events after each cut for the data and the expected background for the 
missing energy channel. The background estimates are normahzed to 10.4 pb~^. The quoted 
error is statistical. The last column shows the selection efficiencies for the uulU^—* all) final 
state, for a 68 GeV Higgs boson. 



5 The Tau Channels 

The r+r^qq final state can be produced via the processes e+e^— s>Z°H°— i>r+r^qq (about 3% of 
the total Z°H° production rate) and qqr+r~ (about 5.6%). This analysis is sensitive to both 
processes, which are characterized by a pair of tau leptons and a pair of energetic hadronic jets. 
In addition, either the pair of hadronic jets or the pair of tau leptons should have an invariant 
mass consistent with the Z° mass. These characteristics are used to suppress the backgrounds, 
predominantly from Z^/'j*—^qq and four-fermion processes. 
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The selection begins with the identification of tau leptons using three algorithms which 
address the different decay channels of the tau lepton. 



(a) An electron, identified by a neural network algorithm is classified as a r^— e^uu 
decay if its momentum is larger than 2 GeV, and it is isolated. In particular, the number of 
electromagnetic clusters within a cone of 26° half-angle around the electron track, N^^, must 
be less than six, and the ratio of the electromagnetic energy within an 11° cone to that within 
a 30° cone, RlK^^ must be greater than 0.7. There must be no hadronic calorimeter cluster 
with energy greater than 0.6 GeV that is associated with the electron track. Electrons from 



photon conversions are rejected using a neural network algorithm p6 |. 



(b) A muon, identified using standard selection algorithms is classified as a r^— >yU^z/z/ 
decay if its momentum is larger than 3 GeV, and if it is isolated. In particular, A^^^^ < 5, and 
the ratio of the scalar sum of all track momenta within an 11° cone to that within a 30° cone, 

must be greater than 0.7. 

(c) The remaining tau lepton decays are identified as narrow, isolated jets. Jets are re- 
constructed using a cone algorithm |^ with a half-angle of 23° and with at least 3 GeV of 



associated energy. Within each resulting jet, the sub-jet of 11° half-angle having the highest 
energy is formed. The sub-jets are accepted as tau candidates if they satisfy the fiducial re- 
quirement I cos^l < 0.92, have one or three associated tracks, have an invariant mass less than 
3.5 GeV, and are isolated with Rl^^^ > 0.6. 

In the selection that follows, the tau lepton momentum is approximated by the momentum 
of the visible decay products. When there are two tau lepton candidates with momentum 
vectors separated by less than 23°, one being identified as a leptonic decay (algorithms (a) or 
(b)) and one as a narrow jet (algorithm (c)), only the candidate identified as a leptonic decay 
is selected. 



(1) Events are required to have at least two tau lepton candidates, each with charge of |g| = 1. 

(2) The total track multiplicity of the event must exceed eight. 

(3) Most of the two-photon and e^e^^Z°7 background events are eliminated by requiring 
that the energy in the forward detector, gamma catcher, and silicon-tungsten luminometer 
be less than 4, 10, and 10 GeV, respectively, that | cos6'miss| < 0.97 and that > 3 GeV. 
In addition, the scalar sum of all track and cluster transverse momenta is required to be 
larger than 40 GeV. 

(4) The remaining Z°/7*— >qq background, with and without photon radiation, is further 
suppressed by requiring that events contain at least four jets, reconstructed using the 
cone algorithm with a 23° half-angle as in (c) above (single electrons and muons from tau 
lepton decays are recognized as low- multiplicity "jets"). Events with an energetic isolated 
photon^ are removed. 

*An energetic isolated photon is defined in this context as an electromagnetic cluster with energy larger than 
15 GeV and no track within a cone of 30° half-angle. 
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(5) In signal events, the algorithms (a), (b), and (c) identify 2.3 tau candidates per event on 
the average. Fake tau candidates containing 3 charged tracks are removed by requiring 
that the tracks originate from a common vertex in three dimensions, as reflected by a 
probability of a vertex fit larger than 1%. Fake candidate pairs are further removed by 
requiring that the sum of the track charges be zero and that the tau candidates satisfy 
a pairwise isolation requirement, | cosai ■ cos 02! < 0.8, where is the angle between 
the direction of the ith tau candidate and that of the nearest track not associated with 
it. In those rare instances where more than one candidate pair passes the selection, we 
give preference to those pairs having taus which were identified through their leptonic 
decay mode and, in case of further ambiguity, we choose the pair with the lowest charged 
multiplicity and highest isolation. 

The hadronic part of the event, obtained by excluding the tracks and clusters from the 
selected tau candidate pair, is then split into two jets using the Durham algorithm. The 
invariant masses of the tau lepton pair, mrr, and of the hadron jets, rrihad, are calculated 
using only the tau lepton and jet momentum directions and requiring energy and momentum 
conservation. At this point the selection separates into two parts, one (A) sensitive to the 
(Z''^r+r^)(II°^qq) process and another (B) sensitive to the (Z°— i>qq)(II°^r+r^) process. 

(6) (6A) The selected events must satisfy 75 GeV < nirr < 105 GeV and mhad >30 GeV. 
In addition, E^is is required to be less than 155 GeV, since the neutrinos from the tau 
lepton decays, originating from the Z° boson, give rise to a relatively large missing en- 
ergy. Finally, cuts are implemented to suppress specific four-fermion backgrounds, from 
e"'"e~— i>Z°/7* + Z°/7* and e+e~— >-Z°e"'"e~. If the tau lepton candidates are both classified 
as r^— i> e^uu or both as r^— /i^i/z/, their opening angle is required to be larger than 90° 
and, in the first case, neither electron is allowed to lie within 36° of the beam axis. 

(6B) The selected events must satisfy 75 GeV < mhad < 105 GeV and > 30 GeV. 
Since in this case the mass cuts are less effective against the background, the requirements 
on the properties of the tau lepton candidates are tightened. The opening angle of the 
tau lepton pair must be larger than 110° and, if one of the tau candidates has a track 
multiphcity exceeding two, the pairwise isolation cut is tightened to | cosai-cosa2| < 0.55. 
Furthermore, to suppress four-fermion backgrounds, pairs with leptons of the same flavor 
are rejected. Finally, to suppress the process W'^W"— i>£z/qq, events are rejected if they 
contain any track or cluster with a momentum or energy exceeding 40 GeV. 

Distributions of | cosai ■ cos 02 1 and m^,- are shown in Figures |^ and ^ for the data, the 
backgrounds, and for a 68 GeV Higgs boson signal, with H°— i>r"'"r~. The numbers of observed 
and expected events after each stage of the selection are given in Table |^ where a good agreement 
between data and MC can be seen. The detection efficiency for a 68 GeV Higgs boson is also 
given. No candidate event is observed while the total background for selections A and B is 
estimated to be 0.59ib0.04(stat)±0.14(syst) events. The background systematics is dominated 
by the error due to the modeling of the cut variables. 
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The detection efficiencies for cases A and B as a function of the Higgs boson mass are 
given in Table ^. and include a small correction coming from accelerator-related backgrounds 
in the forward detectors (2.3 %) which are not fully simulated. The detection efficiencies are 
affected by the following uncertainties: Monte Carlo statistics, 2.8%; uncertainty in the tau 
lepton identification efficiency, 4.3%; uncertainties in the modeling of cut variables excluding 
the tau lepton identification, 9.1% (case A) and 7.6% (case B); uncertainties in the modeling of 
fragmentation and hadronization, 1.2%. Taking these uncertainties as independent and adding 
them in quadrature results in a total systematic uncertainty of 10.5% (case A) and 9.2% (case 
B) (relative errors). 



Cut 


Data 


Total bkg. 


qq(7) 


4-ferm. 


77 




e(%), case A 
68 GeV 


e(%), case B 
68 GeV 


(1) 


857 


611.5 


77.8 


79.9 


421.6 


32.2 


58.5 


59.3 


(2) 


358 


306.9 


75.2 


36.8 


194.9 





58.3 


58.6 


(3) 


50 


55.1 


23.6 


31.2 


0.3 





54.0 


52.9 


(4) 


37 


40.1 


15.2 


24.7 


0.2 





51.7 


50.8 


(5) 


15 


20.1 


6.9 


13.2 


< 0.07 





41.8 


40.8 


(6-A) 





0.41 ±0.03 


< 0.01 


0.41 


< 0.07 





22.9 




(6-B) 





0.18 ±0.02 


0.02 


0.16 


< 0.07 







18.9 



Table 3: The numbers of events after each cut for the data and the expected background for 
the tau channels. The background estimate is normalized to 10.4 pb~^. The quoted errors are 
statistical. The last two columns show the selection efficiencies, for cases A and B, for a 68 GeV 
Higgs boson. 



6 The Electron and Muon Channels 

The i~^i~qq {i = e or fi) final states arise mainly from the process e"'"e~— >-Z°H''— i>£"'"£~qq. They 
amount to approximately 6% of the Higgs boson signal topologies with a small contribution 
(0.1% for mH0=68 GeV) from the Z°Z° fusion process e"'"e~— i>e"'"e~H'^— i>e"'"e~qq. 

The analysis adopted concentrates on those final states proceeding through the ffist process. 
These yield a clean experimental signature in the form of large visible energy, two energetic, 
isolated, oppositely-charged leptons of the same species reconstructing to the Z° boson mass, 
and two energetic hadronic jets. The dominant backgrounds are Z°/7*— >-qq and four-fermion 
processes. 

The selection proceeds as follows: 

(1) The selected events are required to have at least six tracks and are reconstructed as 
four jets, using the Durham algorithm with a cut of 1/34 > 0.001 (single electrons or 
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muons are considered as low- multiplicity "jets"). The events must satisfy the relations 
Ip^iJ < (E^is - 0.5^5) and E^is > O.Q^/s. 

(2) The selected events must contain at least one pair of oppositely charged, same flavor 
leptons (e or /i). Muon candidates are identified using standard algorithms The 
identification of electron candidates, optimized for high energy electrons, starts by se- 
lecting electromagnetic calorimeter clusters having energy larger than 5 GeV and 90% of 
the cluster energy is deposited in four or five calorimeter cells at most, if the polar angle 
of the cluster satisfies | cos 6^1 < 0.75 or |cos^^| > 0.75, respectively. If such a cluster is 
found, tracks are sought within 5.7° of the cluster, having a momentum larger than 2 GeV 
and normalized energy-to-momentum ratios (-E/p)norm > ~5. If more than one track 
satisfies these conditions, the one forming the smallest angle with the cluster direction is 
considered as the electron track. In the region | cos6\ < 0.85, other clusters within 10° of 
the original one are merged, provided that there is no other track within 20° . 

Muon candidates are grouped into pairs without further restrictions. For an electron pair 
to be considered, at least one electron candidate must have —2.5 < (-E/p)norm < 5. 

If more than one pair of leptons of the same flavor is found, the pair with invariant mass 
closest to the boson mass is considered. 

(3) Both leptons in the candidate pair must have an energy larger than 20 GeV with at least 
one of them larger than 30 GeV. The energy of an electron candidate is obtained from 
the associated electromagnetic clusters while for a muon candidate it is approximated by 
the track momentum. 

(4) The rest of the event, obtained by excluding the candidate lepton pair, is reconstructed 
as two jets using the Durham algorithm. An explicit lepton isolation cut is made to 
reject the remaining background from Z°/7*— i>qq(7) with real or fake leptons close to 
the hadronic jets, by requiring that each of the leptons has a transverse momentum, 
calculated with respect to the nearest jet axis, larger than 10 GeV. The background from 
e"'"e~— >-Z*^/7* + Z°/7* is suppressed by requiring that the opening angle of the jet pair be 
larger than 50°. 

(5) The selected events must have a lepton pair with an invariant mass consistent with the Z° 
boson mass. For electrons the invariant mass of the lepton pair must lie between 75 GeV 
and 105 GeV, while for muons it must lie between 60 GeV and 120 GeV. The differing 
mass windows take into account the differing resolutions for electrons and muons. 



Distributions of the lepton energy and of the lepton pair invariant mass are shown in Fig- 



ures ^ and |I0|, for the data, the simulated background, and for a simulated signal having 
mHO=68 GeV. The numbers of observed and expected events after each stage of the selec- 
tion are given in Table ^, together with the detection efficiency for a 68 GeV Higgs boson. 
After all cuts, no event survives in either the electron and muon channels, while in total 
0.14±0.02(stat)±0.06(syst) events are expected. 

^{E/p)noim — [iE/p) — i]/cr where E andp are cluster energies and track momenta, and a the error associated 
to E/p, obtained from the measurement errors of E and p. 
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The detection efficiencies as a function of the Higgs boson mass are given in Table ^ These 
are affected by the following systematic uncertainties: Monte Carlo statistics, 1.0% (electron), 
0.9% (muon); uncertainties in the electron (muon) identification, 0.5% (0.4%); uncertainties 
in the modeling of fragmentation and hadronization, 0.3%; uncertainties in modeling the cut 
variables excluding lepton identification, 0.4%. Taking these uncertainties as independent and 
adding them in quadrature results in a total systematic uncertainty of 1.2% for the electron 
channel and 1.1% for the muon channel (relative errors). 



Cut 


Data 


Total bkg. 


qq(7) 


4-ferm. 


Efficiency mH0=68 GeV 


electron (%) 


muon (%) 


(1) 


313 


335.9 


246.1 


89.8 


88.9 


82.8 


(2) 


71 


75.1 


52.7 


22.4 


75.1 


78.0 


(3) 


2 


2.0 


0.9 


1.1 


70.5 


74.9 


(4) 


1 


0.7 


0.2 


0.6 


67.9 


72.1 


(5) 





0.14±0.02 


0.02 


0.12 


65.3 


70.3 



Table 4: The numbers of events after each cut for the data and the expected background in 
the lepton channels. Background estimates are normalized to the integrated luminosity. The 
quoted error is statistical. The last two columns show the detection efficiencies for the processes 
e"'"e~— >-fe"'"e~ or ^~^fi~) H° for a 68 GeV Higgs boson. 



7 Mass Limit for the Standard Model Higgs boson 



The signal detection efficiencies and the numbers of expected signal events, as a function of the 
Higgs boson mass, are summarized for all search channels in Table ^. 

The following uncertainties affecting the numbers of expected signal events are common to 
all search channels: the uncertainty in the integrated luminosity: 0.6%; the uncertainty in the 
Higgs boson production cross section |]13|, which includes that from the collider energy: 1%; 



and the uncertainty in the Higgs decay branching ratios: 2% |^ |T^. Taking these uncertainties 
as independent and adding them in quadrature results in a systematic error, common to all 
search channels, of 3% (relative). In estimating the number of expected events for an assumed 
Higgs boson mass, these uncertainties are added in quadrature to those affecting the individual 
search channels. 



To derive a new limit on the Higgs boson mass, this search, with one candidate event in 
the four jets channel at 172 GeV (mH0=75.6 ± 3.0 GeV where the error is on a Gaussian ffi 
to the peak of a 75GeV Higgs MC distribution) is combined with earlier OPAL searches at 
a/s K,mio with one candidate in the leptonic channel (mHO=61.2 ± 1.0 GeV) and 161 GeV [^], 
with one candidate in the missing energy channel (mHO=39.3 ±4.9 GeV). Two candidates from 
earlier searches with mno < 25 GeV are not considered further. The expected numbers of Higgs 
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boson events, from this search and the combination of the present and earher OPAL searches, 
are hsted in the last two columns of Table |^. These numbers are affected by total uncertainties 
of less than 10%. 



rriiio 


-ttO 

qqH^ 


-ttO 


+ — ttO 


-ttO 

qqR^ 


e^e R 




y/S =170- 


Grand 


(GeV) 


RO^bb 




RO^qq 


RO^r+r- 






172 GeV 


total 


40.0 


14.6(1.7) 


34.0(1.4) 


28.3(0.2) 


2.4(0.0) 


59.6(0.4) 


67.0(0.4) 


4.1 


106.1 


50.0 


21.7(2.1) 


42.0(1.4) 


28.7(0.1) 


8.3(0.1) 


63.0(0.3) 


71.0(0.4) 


4.4 


34.8 


55.0 


24.7(2.1) 


46.6(1.3) 


28.0(0.1) 


11.6(0.1) 


63.2(0.3) 


70.1(0.3) 


4.3 


19.6 


60.0 


27.0(1.9) 


46.9(1.2) 


26.6(0.1) 


15.1(0.1) 


63.6(0.3) 


69.4(0.3) 


3.8 


12.1 


65.0 


28.2(1.7) 


44.4(0.9) 


24.5(0.1) 


17.9(0.1) 


64.8(0.2) 


69.7(0.2) 


3.2 


6.7 


67.5 


28.4(1.5) 


43.0(0.8) 


23.2(0.1) 


18.8(0.1) 


65.2(0.2) 


70.2(0.2) 


2.8 


4.7 


68.0 


28.4(1.5) 


42.7(0.8) 


22.9(0.1) 


18.9(0.1) 


65.3(0.2) 


70.3(0.2) 


2.8 


4.4 


70.0 


28.2(1.3) 


41.3(0.7) 


21.7(0.1) 


19.3(0.1) 


65.3(0.2) 


70.6(0.2) 


2.4 


3.2 


75.0 


26.6(0.8) 


34.3(0.4) 


18.2(0.0) 


18.3(0.0) 


63.5(0.1) 


71.3(0.1) 


1.5 


1.6 


80.0 


23.3(0.2) 


19.2(0.1) 


14.0(0.0) 


14.1(0.0) 


57.6(0.0) 


70.8(0.0) 


0.4 


0.4 



Table 5: Detection efEciencies (in %) and numbers of expected Higgs boson events (between 
parentheses) at 172 GeV for each search channel separately as a function of the Higgs boson 
mass. The quoted efficiencies are obtained from a fit to values determined at fixed values of 
m^o . The last two columns show the total numbers of expected events in the present search at 
170 and 172 GeV, and the grand total, which also includes the expectations from earlier OPAL 
searches at center-of-mass energies close to the mass and the revised 161 GeV analysis. 



Figure |TT| shows separately the number of expected events for the present search, for previous 
OPAL searches at the Z° peak and ^/s = 161 GeV (where the revised 161 GeV analysis is 
implemented), and for their sum, as a function of the Riggs boson mass. Also shown is the 
95% confidence level upper limit on the number of observed candidate events. A lower limit on 
the Riggs boson mass, of 69.4 GeV, is extracted at the 95% Confidence Level. In deriving this 
limit, the probability that a candidate event with a given observed mass actually originates 



from a Riggs boson of arbitrary mass is calculated following Ref. |29| which takes into account 
the experimental mass resolution including tails in the various search channels. The expected 
background is reduced by the systematic error per channel and then subtracted. It was found 
that the errors on the background estimation have marginal effects on the results, reducing the 
derived limit by 0.1 GeV. The systematic errors are incorporated into the limit according to the 



method prescribed in Ref. [PD|, reducing the derived limit by an additional 0.1 GeV. The effect 
of channels weighting and background subtraction is small. If all channels were assigned equal 
weights, irrespective of the expected rate and background, and no background subtraction had 
performed, the limit would have gone down to only 69.0 GeV. 

Figure |T2| shows the measured confidence level and the expected one (averaged over a large 
number of hypothetical experiments with no signal and candidates spread according to the 
expected background distributions) as a function of the Riggs boson mass. From this Figure it 
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can be seen that the expected limit is at 65.4 GeV and the expected C.L. for the experimental 
limit of 69.4 GeV is 82%. According to MC trial experiments a reasonable value of 11.8% was 
found for the probabihty of inferring a mass limit greater than or equal to 69.4 GeV assuming 
no contribution from a Higgs boson signal. 



8 Summary 

A new search is presented for the Standard Model Higgs boson produced in association with a 
fermion-antifermion pair. The search is based on data collected in 1996 by the OPAL experi- 
ment at center-of-mass energies of 170 and 172 GeV, with an integrated luminosity of 10.4 pb~^. 
The data show no significant excess beyond the background predicted by the Standard Model. 
Combined with earlier OPAL searches at center-of-mass energies in the vicinity of the Z° reso- 
nance and a revised analysis of the 161 GeV data, this search leads to a lower hmit of 69.4 GeV 
for the mass of the Standard Model Higgs boson, at the 95% confidence level. 
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Figure 1: Four jets channel: distribution of input variables 3,5,6,8,9 and 14 (as described in the 
text) for OPAL data at ^^=170 and 172 GeV compared with the Monte Carlo expectations; 
data: points with error bars, simulations (normalized to recorded luminosity): open / shaded 
/ dashed histograms for Z'^/7*— >qq / four-fermion processes / signal (mHO=68 GeV) scaled by 
a factor of 20. All distributions are shown after the preselection (cuts 1-6). 
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Figure 2: Four jets channel: (a) Likelihood distribution for OPAL data at ^/s — 161 GeV 
compared with the Monte Carlo expectation, (b) Likelihood distribution for OPAL data at 
^/s = 170 and 172 GeV compared with the Monte Carlo expectation. All distributions are 
shown after the preselection (cuts 1-6). 




Figure 3: The candidate event in the four jet channel. Top: r0 view with tracks and clusters; 
bottom: zoom in on the vertex region. The jets associated to the Higgs in the kinematic fit 
yielding the highest probability, are labeled as Higgs jets ('H'-jet); the others are labeled as Z° 
jets ('Z'-jet). The error ellipse of the second Higgs jet is too small to be seen in this plot. 
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Figure 4: Four jets channel: likelihood distribution of OPAL events at ^/s — mzo scaled to 
LEP 2 energies compared with the Monte Carlo expectation. 




Figure 5: Missing energy channel; distributions of the missing mass squared after cut (5) 
data: points with error bars, simulations (normalized to recorded luminosity): open / shaded 
/ dashed histograms for Z°/7*^qq / four-fermion processes / signal (mH()=68 GeV) scaled by 
a factor of 50; The arrows indicate the cut position, 76^ < m^jgg < 120^ GeV^. 
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Figure 6: Missing energy channel; scatter plot of S5 vs. the visible mass, after cut (7); dots: 
simulated four-fermion background; open squares: signal for mHO=68 GeV. straight line: cut 
region. The one event (not shown) which survived the cuts, is at a visible mass of 67.7 GeV 
and E5=1.6. 




Figure 7: Tau channels; distributions of the pairwise isolation parameter (see text) after cut 
(3); data: points with error bars, simulations (normalized to recorded luminosity): open / 
shaded / dashed histograms for Z*^/7*-^qq / four-fermion processes / signal (mH0=68 GeV ). 
Arrows indicate domains accepted by the cuts. 
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Figure 8: Tau channels; distributions of the r+r" invariant mass after cut (4); data: points 
with error bars, simulations (normalized to recorded luminosity): open / shaded / dashed 
histograms for Z°/7*^qq / four-fermion processes / signal (mH0=68 GeV ). Arrows indicate 
domains accepted by the cuts. 
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Figure 9: Electron and muon channels; scatter plot of the energies of the lepton candidates, 
after cut (2). The symbols for the data and the simulated signal (mH0=68 GeV, e+e^ and 
channels combined) as indicated; small dots: simulated backgrounds. Solid lines: position of 
cut (3). 
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Figure 10: Electron and muon channels; invariant mass of the lepton pair, after cut (4). 
Open / shaded / dotted histograms: Z°/7*— s>qq / four-fermion backgrounds / simulated signal 
(mH0=68 GeV, e+c^ and fi^fi^ channels combined). All simulations are normalized to the 
integrated luminosity. Horizontal arrows: selected mass ranges; vertical arrow: invariant mass 
of the event remaining in the e"'"e~ channel after cut (4). 
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Figure 11: Expected number of events as a function of the Higgs boson mass. OPAL searches 
at ^/s—mzo (dotted), including y^=161 GeV (dashed) and the grand total (full line) including 
OPAL searches at -v/i=170/172 GeV. The dash-dotted horizontal line is the 95% confidence 
level upper limit for a possible Higgs boson signal in the presence of the observed high-mass 
candidate events (see text). The intersection of the solid curve with the dash-dotted curve, 
indicated by the arrow, determines the 95% confidence level lower limit obtained for the Higgs 
boson mass; 
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Figure 12: The expected (dashed) and measured (sohd) Confidence Level (1-C.L.) as a function 
of the Higgs boson mass when all channels at all center-of-mass energies are combined. The 

intersection of the 95% CL. line (dash-dotted) with the dashed and solid curves determines the 
95% expected and measured (respectively) confidence level lower limits obtained for the Higgs 
boson mass; the latter is indicated by an arrow. 



